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A stability criterion is derived in general relativity for self-similar solutions with a scalar field 
and those with a stiff fluid, which is a perfect fluid with the equation of state P — p. A wide 
class of self-similar solutions turn out to be unstable against kink mode perturbation. According 
to the criterion, the Evans-Coleman stiff-fluid solution is unstable and cannot be a critical solution 
for the spherical collapse of a stiff fluid if we allow sufficiently small discontinuity in the density 
gradient field in the initial data sets. The self-similar scalar-field solution, which was recently found 
numerically by Brady et al. (2002 Class. Quantum. Grav. 19 6359), is also unstable. Both the flat 
^— . Friedmann universe with a scalar field and that with a stiff fluid suffer from kink instability at the 

^— ■>) ' particle horizon scale. 
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I. INTRODUCTION 

Since the Einstein equations are simultaneous nonlinear partial differential equations, it is not easy to obtain general 
solutions. In this context, we can assume self- similarity and consider self-similar solutions. For example, in spherical 
symmetry, this assumption reduces the partial differential equations to a set of ordinary differential equations and 
even the classification of self-similar solutions with a perfect fluid has been done H, Q based on the theory of 
dynamical systems ([^|^ and references therein). Thus, the assumption of self-similarity is very powerful in finding 
dynamical and inhomogeneous solutions. The application of self-similar solutions is very large, including cosniological 
perturbati ons, star formation, gravitational collapse, primordial black holes 0, S j cosmological voids and cosmic 
• censorship [lOl [ill . See ^3 for a brief review of self-similar solutions. 

\ In recent progress in general relativity, self-similar solutions have attracted attention not only because they are easy 
O . to obtain but also because they may play important roles in cosmological situations and/or gravitational collapse. One 
^ '■ of the important roles of self-similar solutions is to describe asymptotic behaviours of more general solutions. This 
$H feature has been found in several models such as homogeneous cosmological models [^|^. In a cosmological context, 
the suggestion that spherically symmetric fluctuations might naturally evolve from complex initial conditions to a 
^ , self-similar form has been termed the similarity hypothesis |l3l Il4| . The hypothesis has been also suggested to hold 
! ■ in more general situations, including gravitational collapse. In fact, Harada and Maeda numerically found that a 
rS ' self-similar solution acts as an attractor in the spherically s ymm etric collapse of a perfect fluid with the equation of 
^ . state F = /cp (0 < /c < 1) at least for < fc < 0.03 (see also"|ll|) 
■ ■ ' Spherically symmetric self-similar solutions have been also studied in the context of critical behaviour in gravitational 
collapse, in which a self-similar solution is not an attractor but an intermediate attractor. The critical behaviour was 
numerically discovered by Choptuik 17] in the spherical system of a massless scalar field. There appears a discretely 
self-similar solution at the threshold of black hole formation for p = p* , where p and p* are the parameter of initial 
data sets and the critical value of the parameter, respectively. For supercritical collapse, the black hole mass AfsH 
follows the scaling law Mbh cx: \p — p*\^ , where 7 « 0.37 is referred to as the critical exponent. Subsequently, these 
findings have been confirmed by several authors (e.g. 10). Evans and Coleman studied the spherical system of a 
radiation fluid (P = p/3) and found similar phenomena although the observed critical solution is a continuously self- 
similar solution. Koike, Hara and Adachi j2y| showed that the critical behaviour is understood by the intermediate 
behaviour around a saddle equilibrium point. This approach was also successfully applied to the spherical system of a 
massless scalar field to understand the critical behaviour observed by Choptuik 21]. Several authors i15b,2^.22] have 
extended the work to perfect fluids with the equation of state P = fcp (0 < fc < 1). The limiting case fc ^ turns out 
to be the Newtonian critical phenomena Another case fc = 1, for which the perfect fluid is called a stiff fluid, is 
interesting because the stiff fluid is known to be equivalent to a massless scalar field in certain circumstances (2^ .26j . 
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Neilsen and Choptuik '23l| studied the stiff-fluid system and found that the critical behaviour can be regarded as a 
continuous Hmit of that for a perfect fluid with P = fcp (0 < A: < 1) as fc ^ 1. The critical solution is continuously 
self-similar and the critical exponent is given by 7 = 0.92 ± 0.02. Brady et al. l2fl studied the mystery of why 
the critical behaviour for a stiff fluid is different from that for a scalar field. See [23| for a broad review of critical 
phenomena in gravitational collapse. 

The present article concerns weak discontinuity around self-similar solutions. It has been known that the required 
differentiability condition for the spacetime metric tensor can be much weaker than C°° , for the existence and unique- 
ness of geodesies and/or well-defined initial value formulation of the Einstein equation Moreover, it has been well 
established that shock waves, which is one of the examples of strong discontinuity, play essential roles in many astro- 
physical situations. If we insert a perturbation with a sufficiently small weak discontinuity, it may grow too large to 
neglect as time proceeds. This instability has been called kink instability. The kink instability of self-similar solutions 
was studied by Ori and Piran 29] for the spherical system of isothermal gas in Newtonian gravity. Harada (paper 
I) investigated the kink instability of self-similar solutions for the spherical system of a perfect fluid with the equation 
of state P = fcp (0 < fc < 1) in general relativity. In this article we derive a stability criterion against the kink mode 
perturbation for both self-similar solutions with a scalar field and those with a stiff fluid in general relativity. Several 
interesting applications are discussed, one of which is related to critical behaviour for these two systems. 

This article is organized as follows. In Sec. II, we derive the basic equations for the spherically symmetric Einstein- 
scalar-field system and those for self-similar solutions. Also in this section we introduce the kink mode for the 
self-similar scalar-field solutions, derive the equations for this mode and obtain the stability criterion for self-similar 
scalar-field solutions. In Sec. Ill, the analysis in paper I is extended to include the stiff-fluid system. Section IV is 
devoted to obtain the correspondence relation between self-similar solutions of the two systems, the scalar-field system 
and the stiff-fluid one. In Sec. V, we discuss the applications of the present result, in particular, to critical behaviour 
of gravitational collapse. In Sec. VI, we summarize the paper. We adopt units such that G = c = 1 and the abstract 
index notation of 31,]. 

II. STABILITY CRITERION FOR SELF-SIMILAR SOLUTIONS WITH A SCALAR FIELD 
A. Equations for a self-gravitating massless scalar field 

In this section, we basically follow the notation given by Brady et al. [2^. We consider a massless scalar field as 
a matter field, for which the stress-energy tensor is given by 

T'^'' = V'^./.V''./.- i<7"''(Ve0V'=0). (2.1) 

We adopt the Bondi coordinates for a spherically symmetric spacetime as 

ds^ = -ggdu^ - 2gdudR + R^idO^ + sin^ Odcj)^), (2.2) 

where g = g{u, R) and g = g(u, R). Then the Einstein equation and the equation of motion for the scalar field reduce 
to the following partial differential equations: 

i\ng),R = AnR^%, (2.3) 
R9,R ^ g-9, (2.4) 

g(^^^ = 87rP((/)2„-g0,„0,^), (2.5) 

{R''gcj)M),n. = 2R^^u + 2R^(j)^uR, (2.6) 
where the comma denotes the partial derivative. For later convenience we define a new function h{u, R) as 

(p ^h{u,R) - K\n\u\, (2.7) 
where k is an arbitrary constant, and we define the following variables: 

x = ~-, X = \n\x\, T=-\n\u\. (2.8) 
u 
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Then equations H2.3|l - H2.6|l reduce to 



(In.g)' = 47r/i'2, (2.9) 
g' = 9-9, (2.10) 



9_ 

9 J \9 



= 8TT{h + h' + K)[x{h + h' + k) ~ gh'], (2.11) 
{fh'y + fh' = 2x[{h + h') + {h + h')' + k], (2.12) 



where the dot and prime denote the partial derivatives with respect to T and X, respectively. We refer to it < 
and M > as early times and late times, respectively. We refer to u < and u > as early times and late times, 
respectively. It is noted that the present definition of prime is different from that in Brady et al. [2^. 

B. Self-similar scalar- field solutions 

For self-similar solutions, we assume that g — g{X), g — g{X) and h = h{X). The ordinary differential equations 
for g, g and h are given by 

{\ng)' = 4^/, (2.13) 

9' = 9-9, (2.14) 

g-g ^ 47r[2^x-{g~2x){f + 2Kj)], (2.15) 

h' = j, (2.16) 

{2x-g)j' = -2>ix + jig-2x), (2.17) 

where the prime denotes the differentiation with respect to X . 

These equations are singular when g — 2x, which is called a similarity horizon. We denote the value of x at the 
similarity horizon as Xs and also the value of X as Xs = In |xs|. This line corresponds to a radial ingoing (outgoing) 
null curve for early-time (late-time) solutions because 

along a radial ingoing (outgoing) null curve. We consider self-similar scalar-field solutions with finite values of functions 
g, g, h and j and their gradients with respect to X at the similarity horizon. Then we find sd X — X^ 

g = 2xs{ATr^ + l), g = 2xs, J = t^, (2.19) 

inn 

for K and 

g = g = 2x3, j = is, (2.20) 
for K — 0. For k = 0, similarity horizons are parametrized by js- Since we have 



2x) 



1-4W, (2.21) 



at the similarity horizon, the self-similarity horizon is "transluminal" , i.e. it has subluminal interior and superluminal 
exterior for Attk^ < 1, while it is "anti-transluminal" , i.e. it has superluminal interior and subluminal exterior for 
Attk"^ > 1. 

We can have regular centre at i? = for u ^ 0, which corresponds to X = — oo. The regularity condition reduces 
to the following initial condition at X = — cx). 

g = g = go, h = ho, J = 0- (2.22) 
Taking the limit X — oo in Eq. (|2.17|) and using the I'Hospital's theorem, we have 

x^-oo x{X) x^-oo x{X) go 
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The above system of Eqs. (|2.13|) " (|2.17|) was studied b y C hristodoulou [33 |. Brady 'ss'l and Brady et al. Exact 
self-similar scalar-field solutions discovered by Roberts [m| correspond to general solutions for k = 0, which will be 
described in some detail in Appendix IXI Christodoulou found that there is a solution with 47rK^ = 1/3, which has 
both regular centre and analytic similarity horizon and turns out to be the flat Friedmann solution. Using a two-point 
shooting method, Brady et al. 26] numerically found another self-similar scalar-field solution with both regular centre 
and analytic similarity horizon, which we will discuss in Sees. IV^ and IVBl in detail. 



C. Kink instability of self-similar solutions with a scalar field 

We consider perturbations which satisfy the following conditions in the background of a self-similar solution: (1) 
The initial perturbations vanish inside the similarity horizon for early-time solutions (u < 0). (For late-time solutions 
(u > 0), the initial perturbations vanish outside the similarity horizon.) (2) g, g, h and h' are continuous everywhere, 

in particular at the similarity horizon. (3) h" and h" are discontinuous at the similarity horizon, although they have 
finite one-sided limit values as X — > — and X ^ + 0. 
We denote the full-order perturbations as 

Sg{T,X) = g{T,X)-gi,{X), 5giT, X) = g{T, X) - gb{X), 5h{T , X) = h{T , X) - h^{X) , (2.24) 

where g^,, and /ib denote the background self-similar solution. 

By conditions (2), the perturbations satisfy 5h = Sh' = and Sh" ^ at the similarity horizon at initial moment 
T = Tq. The evolution of the initially unperturbed region is completely described by the background self-similar 
solution because no information from the perturbed side can penetrate the unperturbed side by condition (1). Then, 
by conditions (2) and (3), we find Sh = Sh' = and Sh" 7^ at the similarity horizon for T > To for early-time 
solutions (for T <Tq for late-time solutions). We find 5g — Sg = Sg' = Sg' — 5g" ~ but Sg" 7^ and Sg'" ^ from 
Eqs. (|2.9|) . H2.10|l and (|2.11|) for the full-order perturbations at the point of discontinuity. 

Differentiating Eq. H2.9|l with respect to X and estimating both sides at the point of discontinuity, we obtain 

Sg" = 87r5bJb<5/. (2.25) 

Using Eqs. 12.9|l and H2.10() . differentiating Eq. (|2.11() with respect to X, and estimating both sides at the point of 
discontinuity, we obtain 

(jb + «)(5b-2a;)5/ = 0. (2.26) 

Therefore, Sj' 7^ is possible at the similarity horizon X = Xs where g^, — 2xs. 

Using Eqs. H2.9|l . H2.1()|l and H2.26() . differentiating Eq. H2.11|l twice with respect to X , and estimating both sides at 
the point of discontinuity, we obtain 

2a;(jb + n)Sj' ~ [ffbib + 2(.gb - 4a;)(jb + >i)]Si' = (jb + «)(.% - 2x)5j" . (2.27) 

We can have another equation by using Eq. H2.1(J|I . differentiating Eq. H2.12|l with respect to X, and estimating both 
sides at the point of discontinuity as 

2x5]' - (2gb - 5b - 4.x)6j' = (5b - 2x)Sj". (2.28) 

The two equations obtained above are independent in general but degenerate at the similarity horizon. At the 
similarity horizon, we find that these two equations reduce to 



Sj' - (87r^2 _ ^ ^li^^ _ 1^ ^^..] (2.29) 



where Eq. H2.19|l or Eq. (|2.2U|) is used. Since the left-hand side has a finite limit at the similarity horizon, the right- 
hand side also has a finite limit. Since Sj" may not have a finite limit at the similarity horizon, the limit value on 
the right-hand side is not trivial. However, we can show that the right-hand side vanishes employing the proof in 
Appendix Therefore, we obtain the full-order perturbation equation for Sj' at the similarity horizon: 

Sj' - (Bttk^ - l)Sj' = 0. (2.30) 

It should be noticed that the full-order analysis results in a linear-order equation. This equation can be integrated as 

Sj' = const • e"'^, (2.31) 
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where 

a = -(1 - SvTK^). (2.32) 

Therefore, for early-thiie solutions, it is found that the perturbation decays exponentially for Attk^ < 1/2, it is constant 
for Attk^ = 1/2 and it grows exponentially for Attk'^ > 1/2. The situation is reversed for late-time solutions. 

Here we define instability by the exponential growth of discontinuity. Then we find the following criterion: for 
early-time solutions, solutions with a regular similarity horizon and inn^ < 1/2 are stable against the kink mode, 
while those with Attk'^ > 1/2 are unstable. Solutions with 4ttk'^ = 1/2 are marginally stable against this mode. The 
situation is reversed for late-time solutions. 



III. STABILITY CRITERION FOR SELF-SIMILAR SOLUTIONS WITH A STIFF FLUID 



A. Equations for a self-gravitating stiff fluid 

In Sec. Ill Cl we have seen that the kink instability can occur in self-similar scalar-field solutions. Since a scalar 
field is equivalent to a stiff fluid, it is expected that self-similar stiff-fluid solutions can also have the kink instability. 
However, the equivalence turns out to be incomplete as we will see in Sec. IIVI Hence it will be interesting to see 
when and how the kink instability occurs in stiff-fluid self-similar solutions. In paper I, the kink instability is studied 
in perfect-fluid systems with the equation of state P = kp {0 < k < 1). Almost all equations in paper I are available 
also for a stiff fluid with k — 1. Here we avoid repetition of the derivation and only quote the equations which are 
necessary for the subsequent analysis. We basically follow the formulation given by Q. 

The stress-energy tensor for a stiff fluid is given by 

T'''' = p(2u''u'' + .g°''), (3.1) 

where p is the energy density and m° is the four velocity of the fluid element. In a spherically symmetric spacetime, 
there is a natural quasilocal mass called the Misner-Sharp mass, defined by |35j 

m=^il- V'^i?V„i?) , (3.2) 

where R is the circumferential radius. The line element in a spherically symmetric spacetime is given by 

ds^ = -e^'-''''^ dt^ + e'^'-^'^'^dr^ + R^{t,r){d9^ +sm^ 9d^''), (3.3) 

where we choose the radial coordinate r to be comoving with fluid elements. We define dimensionless functions such 
as 

77 = 87rrV, S=-, M=— , y=^, (3.4) 
r r rib" 

where we have also defined an auxiliary function y for later convenience. We introduce the following variables: 

^=^, r = -ln|t|, z = -ln|e|. (3.5) 

The velocity function V is defined by the velocity of the ^ = const curve relative to the fluid element, which is written 
as 

V = sign(Oe'^+'^. (3.6) 
We can find that the equations of motion for the fluid are integrated as 

= a,77-le2^ = a^r^S'^, (3.7) 
where we make Oo- and a^j constant using coordinate transformation. Using Eqs. (|3.7|l . V"^ is rewritten as 

y2 ^ '^s-\ (3.8) 
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Using the dimensionless variables, we can derive the following equations from the Einstein equations: 

M + M' = T]S^{S + S'), 
M + M' = -7]S^{S + S'), 

^ \s + S') 



S + 2S' + S" + 




('1, 


---') 


+ 1 








V J 





V 



-2\{S + S')--e 



52 



and 



M 



-cr+22 / 



(3.9) 
(3.10) 



(3.11) 



(3.12) 



where e"^ and are given by Eq. (|3.7|l . the dot and prime denote the partial derivatives with respect to r and z, 
respectively, and three of the above four equations are independent. 



B. Self-similar stiff- fluid solutions 



For self-similar solutions, we assume that all dimensionless quantities depend only on 0, i.e., M = Af(z), S — S{z), 
rj = r\{z), a — cr{z), and cu = uj{z). Then, from equations H3.9|l - H3.12|l . we obtain the following ordinary differential 
equations: 



M' 
S' 

v' 



I I -y 
2 y 



M, 
5, 



y - agV^ 

1- y2 



and constraint equation 



M 



0, 



(3.13) 
(3.14) 
(3.15) 

(3.16) 



where V'^ is given by Eq. H3.8|) . It should be noted that does not change under the scale transformation and 
therefore its value has a physical meaning. Oo- is physically characterized by the following relation: 



A sonic point z = z^is defined by 



1. 



(3.17) 



(3.18) 



Since the sound speed is equal to the light speed, the sonic point may correspond to the Cauchy horizon, event horizon 
or particle horizon, simultaneously. From Eqs. H3.8|l . (j3.16|l and H3.18|l . we find M = Ms and S = Ss B.t the sonic 
point, where 



1 + flcr \a. 



1/4 



Ss 



1/4 



(3.19) 



Here we adopt the finiteness of the density gradient as a regularity condition. Then, from Eq. H3.15|l . we obtain y = ys 
at the sonic point, where 



ys = a<T- 

We introduce a new independent variable u, which is defined as 

dz o 
— = 1-V^. 
du 



(3.20) 



(3.21) 
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Using u in place of z, we have the following system of the ordinary differential equations: 

'-{l-V'^)M, (3.22) 



du 2 y 

f ^ -i^(l-V-)5, (3^23) 
^ = 2 [(1 - j)(l - F^) + {!, - «,V')] „. (3.24) 

In the above, the sonic point is found to be an equilibrium point of the dynamical system with four dimensional 
state space {(M, S", 77, 2)}. More precisely, the whole of these sonic points composes the two dimensional surface of 
equilibrium points determined by Eqs. (|3.18|l and (|3.2U|) . which is called a sonic surface. The qualitative properties 
of solutions are analyzed by the linearisation of the ordinary differential equations at the equilibrium point. We put 

M = risSl{ys + xi), S = Ss{l + X2), 7] ^ r]s{l + X3), S, ^ (^^{l + X4) (or z^Zs-Xi), (3.25) 

where Xi to X4 are regarded as components of the vector x. Then, we find the following linearised ordinary differential 
equations in the matrix form: 

-^x^Ax, (3.26) 



du 



where the matrix A is given by 




(3.27) 



This matrix has two zero eigenvalues and two generically nonzero eigenvalues. The former two zero eigenvalues are 
due to the two dimensional surface of equilibrium points in four dimensional state space. The latter two eigenvalues 
are: 

Ai=-2y„ and A2 = -2(1 - j/,), (3.28) 
associated with the following eigenvectors: 




ei = -, , and 62 



/ (l-2y,)(l-2/3) 
-(l-2y3)(l-ys) 

V -2(1 -2y,) 



(3.29) 



respectively. Therefore, as we restrict our attention to the plane spanned by ei and 62, the sonic point is a node 
with primary direction ei and secondary direction 62 for < < 1/2, a degenerate node for ys — 1/2, a node with 
primary direction 62 and secondary direction ei for 1/2 < j/g < 1, and a saddle with repulsive-eigenvalue direction 62 
and attractive-eigenvalue direction ei for 1 < j/g- 

Both eigenvectors are regular in the whole four dimensional state space {(M, S, 77, z)} but not in the physical space. 
Actually, along eigenvector ei, the density gradient takes the following limit: 

v' 

lim — = ±00. (3.30) 

z — >Zs Tj 

The above divergence directly implies the divergence of the density gradient with respect to the physical length I on 
the t = const hypersurface as seen in the following equation: 



dl ) t ydrSirr^ J^ Sirr^ \ t] 



(3.31) 



Hereafter, we restrict our attention to the self-similar solutions with finite density gradients at the sonic point. Under 
this restriction it is found that only the solutions along 62 are allowed and that there is no acceptable solution with 
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Us = 1/2. This is quite different from the case for < fc < 1, in which both directions associated with the two 
eigenvalues are physically acceptable. Along direction 82, we have 

f (3.32) 

at the sonic point for 7^ 1/2. Then, from Eq. (|3.8I) . we have 

(V^y = 2il-y,), (3.33) 

at the sonic point. Therefore, the sonic point is "transsonic" for Q < < 1/2 or 1/2 < ys < I, while the sonic point 
is "anti-transsonic" for 1 < y^. 

The classification is as follows. For < j/s < 1/2, the solution crosses a nodal sonic point along a secondary direction. 
For ys = 1/2, the solution crosses a degenerate-nodal sonic point. For 1/2 < j/g < 1, the solution crosses a nodal sonic 
point along a primary direction. For 1 < j/g, the solution crosses a saddle sonic point along a repulsive-eigenvalue 
direction. As for the family of acceptable solutions, for < ?/s < 1/2 or 1 < ?/s there is only one acceptable solution 
which crosses the sonic point and it has analyticity there (see paper I for the analyticity of self-similar solutions at 
the sonic point, [s^) For 1/2 < ?/s < 1, there is a one-parameter family of acceptable solutions which crosses the sonic 
point and only one of them has analyticity there. For ys — 1/2, there is no acceptable solution. It is clear that the case 
for ys = I needs more careful treatment because the equilibrium point is nonhyperbolic (c.f. the Hartman-Grobman 
theorem). Bicknell and Henriksen "o] found that there is a one-parameter family of regular solutions in this case by 
a direct integration around the sonic point. 



C. Kink instability of self-similar solutions with a stiff fluid 

We consider perturbations which satisfy the following conditions in the background of a self-similar solution: (1) 
The initial perturbations vanish inside the sonic point for early-time solutions {t < 0). (The initial perturbations 
vanish outside the sonic point for late-time solutions {t > 0).) (2) M, S and rj are continuous everywhere, in particular 
at the sonic point. (3) 77' and rj' are discontinuous at the sonic point, although they have finite one-sided limit values 
as z — > Zg — and z ^ Zs + [33 . 

We denote the full-order perturbations as 

SS{t, z) = S{t, z) - S^{z), SM(t, z) EE M(t, z) - Mb(z), <5?7(t, z) = rj{T, z) - rn,(z), (3.34) 

where Afb, Sh and 77b denote the background self-similar solution. 

For the perfect fiuid case for the equation of state P = kp, the evolution equation for the discontinuity at the sonic 
point becomes (paper I but 38]) 



77b V l + ^»?b i + k J rib 1 + k \ rjh 



1 f± \ Srf 

2 J 7?b 



(3.35) 



where the finite limit value on the right-hand side exists because the left-hand side has a finite limit value. Employing 
the proof given in Appendix we can show that the limit value on the right-hand side is zero. Then, we have 
(Eq.(4.15) of paper I) 

W I ( qI-M; I 5-3k-4ys Yv' l-k/5v'y _^ 
Vh \ l + kT]b 1 + k J T]b l + k\ r]b J 

With k — 1, we find that the equation for full-order perturbation at the sonic point is given by 

^ + (1-2,^)^=0. (3.37) 

It should be noticed that the full-order analysis results in a linear-order equation. This equation is integrated as 

^ = const • (3.38) 
Vh 

where 



a = -{l-2ys). 



(3.39) 
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Therefore, it is found that as r increases, the discontinuity in 77' decays exponentially for < j/s < 1/2 but grows 
exponentially for 1/2 < yg. We should note that the parameter j/g can be identified with by the discussion given 
in Sec. HTlBl 

This result for the stiff- fluid case is quite different from that for < fc < 1 in the following respect. For the latter 
case, the instability results in the divergence of Srj' for a finite value of r because of its nonlinear growth. However, 
for the former stiff-fluid case, the instability is much weaker and only exponential with respect to t. This weakness of 
kink instability for the stiff-fluid case may make the physical interpretation of instability rather subtle as is mentioned 
in paper I. However, here we define instability by the exponential growth of discontinuity. 

Then we find the following criterion: for early-time solutions, solutions with a regular sonic point and < Us < 1/2 
are stable against the kink mode, while those with a regular sonic point and 1/2 < j/g are unstable against the 
kink mode. The situation is reversed for late-time solutions. If the criterion for kink instability should be related 
to the classification of sonic points, we can give the following criterion: as for early-time solutions, nodal-point 
secondary-direction solutions are stable, while nodal-point primary-direction solutions and saddle-point, repulsive- 
eigenvalue-direction solutions are unstable. The situation is reversed for late-time solutions. The above indicates that 
the stability criterion in terms of the properties of sonic points, which is obtained in paper I for < k < 1, also applies 
to the stiff-fluid case fc = 1. 



IV. CORRESPONDENCE BETWEEN SCALAR-FIELD AND STIFF-FLUID SOLUTIONS 



A. Correspondence between a scalar field and a stiff fluid 



We briefly review the correspondence between a scalar field and a stiff fluid. The stress-energy tensor for the stiff 
ffuid is written by Eq. We assume p > 0. Introducing a scalar field (j) with a future-pointing timelike gradient 

such as 

Va0 = ^/%Ma, (4.1) 

which is possible only for an irrotational velocity field, we have 

p = -iVe(/.V^0, (4.2) 

and it turns out that the stress-energy tensor is written by Eq. (|2.1(l . i.e. that of a massless scalar field. Therefore, 
the irrotational stiff fluid with positive energy density has an equivalent scalar fleld. Actually, the correspondence is 
not one to one but one stiff-fluid solution to two scalar-field solutions because (j) — —4> \s also an equivalent solution 
changing the relation (|4.1|l as 

-Va0 = ^/2'pua. (4.4) 
Inversely, the scalar-field solution has an equivalent physical stiff fluid only if it has a timelike gradient. 



B. Correspondence between self-similar scalar-field solutions and self-similar stiff- fluid solutions 

Here we obtain the correspondence between the self-similar scalar-field solutions studied in Sec. Ill Hi in the Bondi 
coordinates and the self-similar stiff- fluid solutions studied in Sec IIIIBl in the comoving coordinates. 
In the Bondi coordinates (|2.2|l . Eq. 14.2|l yields 

P = (20,« - #,fl) = ^^[22;(i + ^) - 9i]- (4.5) 

The relation between the time coordinates t and u is given by the two equivalent expressions for the proper time 
element at the regular centre as 



(4.6) 
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where the arguments of the left-hand and right-hand sides are z and X, respectively. Integrating the above equation, 
we have 

e'"/2(-c»)t = y/^{-(X))u. (4.7) 
Then we can write the parameter Og- of stiff-fluid solutions in terms of the quantities of scalar-field solutions as 

flo- = lim rye'^e"^^, 



lim 87rp(t, r)r^e'^ 
1 — >o \ r 



2 



= |un 8ttp{u, R)R^gg ^ 
= lim 4TTX~'^gj[2x{j + k) - gj]. 

Therefore, the relation between the self-similar stiff-fluid solution and scalar-field solution is given by 

Qa- = 47rK^ (4.8) 

where Eqs. (|2.22|) and H2.23|l are used. 

Moreover, we can explicitly show that the stiff-fluid solution which is equivalent to the scalar-field solution with 
regular centre and regular similarity horizon satisfies the regularity condition at the sonic point. The definition of the 
Misner-Sharp mass yields 

m^^^fl-l) . (4.9) 



2 V 9 

Then the function y is written in terms of the quantities associated with the scalar-field solution as 

m Q — Q , ^ 

V = — f4 10) 

AirpR^ 4TTj[2x{j + k) - gj]' ^ ' ' 

Then we obtain 

Vs = W, (4.11) 

at the sonic point, where Eqs. H2.19|l or H2.20|l are used. Therefore, the regularity condition H3.20|l at the sonic point 
for stiff-fluid solutions is satisfied. 

Equation (|4.11|) is a crucial correspondence relation between a self-similar scalar-field solution and self-similar 
stiff-fluid solution. As expected, when there is the equivalence between a self-similar scalar-field solution and a self- 
similar stiff-fluid solution, the stability against the kink mode coincides for both self-similar solutions, as is seen from 
Eqs. H2.3HI . 12.32|l . (|3.38|l and l|3.39(l . The growth rates of the kink mode are also the same. 

V. APPLICATIONS 
A. Evans-Coleman stiff-fluid solution and critical behaviour 

There are two important sequences of numerical self-similar solutions with analyticity for a perfect fluid with P — kp, 
the one is the Larson-Penston solution and the other is the Evans-Coleman solution. The Larson-Penston solution in 
general relativity was first discovered and discussed in the context of cosmic censorship by Ori and Piran 10], and 
further analysis was done [Til IT^ l30l|. The solution crosses a node along a secondary direction for < A: ^ 0.036, a 
degenerate node for k « 0.036, and a node along a primary direction for 0.036 ^ k < 1/3. However, the existence of 
the sequence of Larson-Penston solutions for a stiff fluid has not been demonstrated yet. Actually, in our preliminary 
numerical survey, the existence of the Larson-Penston solution was confirmed only for < A; < 1/3. 

The Evans-Coleman solution was first discovered by Ori and Piran identified with a critical solution by Evans 
and Coleman for k = 1/3, and shown to exist and identified with a critical solution for < /c < 1 by several 
authors |23,|U|23. Carr et al. indicated that the character of the sonic point for the Evans-Coleman solution 
changes as follows: the solution crosses a saddle along an attractive-eigenvalue direction for < fc ^ 0.41, a node 
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along a secondary direction for 0.41 ^ fc < 0.89, a degenerate node for k « 0.89, and a node along a primary direction 
for 0.89 < fc < 1. 

For the scalar-field system, Brady et al. numerically solved the set of ordinary differential equations H2.13(l - 
(|2.17l) and searched the value of k for which the solution has regular centre and analytic similarity horizon. They 
found that there is such a solution with Attk^ k, 0.577. We refer to this self-similar solution as the Brady-Choptuik- 
Gundlach-Neilsen (BCGN) solution. They found that this solution has a single unstable mode which is analytic at 
the similarity horizon, and that no physical stiff fluid can reproduce this solution because the BCGN solution has 
a spacelike hypersurface X = Xi{> X^), on which the gradient of the scalar field is null and beyond which it is 
spacelike. They mentioned that it is possible to consider another continuation beyond X = Xi and then the solution 
can be regarded as a physical stiff-fluid solution. Although this solution has some discontinuity ai X = Xi, ii can 
be considered as a continuous limit of the Evans-Coleman perfect-fluid solution with P = kp for 0.89 ^ k < 1 (see 
Fig. 5 of [2^). Hence we refer to this solution as the Evans-Coleman stiff-fluid solution. Neilsen and Choptuik p3|'s 
numerical simulation suggests that this solution would be a critical solution of stiff-fluid gravitational collapse at the 
threshold of black hole formation (see Fig. 2 of |2^). 

First we consider the stability of the Evans-Coleman stiff-fluid solution. The parameter for this solution is given 
by aa = Vs ~ 0.577 using the correspondence relations (|4.8|l and H4.11|) . Then, based on the analysis in Sec. IIIIBl 
the sonic point of the Evans-Coleman stiff-fluid solution is a node and the solution crosses along a primary direction. 
This is quite reasonable for a /c = 1 Evans-Coleman solution since the Evans-Coleman solution crosses a node along a 
primary direction for fc — *■ 1. From the stability criterion obtained in Sec. lIIlCl the Evans-Coleman stiff- fluid solution 
is unstable against the kink mode because ys — 0.577 > 1/2. Actually, the stability against the kink mode does not 
change whatever continuation one chooses beyond X = Xi > X^- The Evans-Coleman stiff- fluid solution has two 
unstable modes, one is analytic and the other is kink. This implies that this solution cannot be a critical solution 
of stiff-fluid gravitational collapse if we allow initial data with sufficiently small discontinuity in the density gradient 
field. It is discussed in paper I why the kink instability did not affect Neilsen and Choptuik |2^'s fiuid simulation. 

B. BCGN scalar-field solution and critical behaviour 

Based on the analysis in Sec. Ill CI the BCGN scalar-field solution is unstable against the kink mode because 
Attk^ ~ 0.577 > 1/2. The BCGN scalar-field solution has two unstable modes, one is analytic and the other is kink. 
This implies that this solution cannot be a critical solution of scalar-field gravitational collapse if we allow initial data 
with sufficiently small discontinuity in the second-order derivative of the scalar field. 

Since this solution was found to have a single analytic unstable mode, this solution is expected to act as an 
intermediate attractor at least locally for some critical behaviour related to the asymptotic behaviour of gravitational 
collapse from the point of view proposed by if we restrict our attention to evolution with C°° functions for 
instance. However, no numerical simulation has been reported so far, that suggests the existence of critical behaviour 
associated with this solution. Brady et al. [2^ stated that the BCGN scalar-field solution is disqualified as a critical 
solution at the threshold of black hole formation because it has an apparent horizon beyond the similarity horizon, 
bmce Choptuik [13 's numerical simulation breaks down at the formation of an apparent horizon because of the 
choice of coordinate system, it would be reasonable that the critical behaviour associated with the BCGN scalar-field 
solution was not observed in his calculation. However, Hamade and Stewart p"s] did not mention the observation of 
critical behaviour associated with a continuously self-similar solution like the BCGN scalar-field solution, although 
their numerical simulation is based on the double-null formulation, which can treat an apparent horizon without any 
coordinate singularity. 

Here we would like to speculate about this puzzle, discrete self-similarity or continuous self-similarity. Of course, we 
should note the possibility that the numerical survey may not have been sufficiently complete yet. If this possibility 
is true, further numerical investigation based on the double-null formulation would reveal a hidden critical behaviour 
associated with the continuously self-similar solution, i.e. the BCGN scalar-field solution, for which the critical 
exponent should be given by 7 = 0.92 ± 0.02. This critical behaviour associated with continuous self-similarity would 
be seen not at the black hole threshold but inside the apparent horizon. In other words, it would be hidden by an 
apparent horizon and therefore by an event horizon and could not be observed by a distant observer. We would have 
two distinct critical behaviours for the same system, one is associated with discrete self-similarity and the other is 
associated with continuous self-similarity. However, we can suggest another possibility that the kink instability, which 
we have analysed here, may prevent the BCGN solution from involving critical behaviour. The possibility depends 
on which class of initial data sets we regard as physically realistic. Mathematically one can avoid the kink instability 
by restricting oneself to solutions in C°° class. In numerical simulations, however, the differentiability condition may 
cause some subtle problems. The above speculations can be proved or disproved by further numerical investigations. 
In each case, it is inferred that the qualitative properties of dynamical solutions in the spherical system of a massless 
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scalar field would be much more complicated than usually supposed in the context of critical behavior. 

C. Exact self-similar solutions 

For the flat Friedmann solution, we find Attk"^ = 1/3 as a scalar-field solution and acr = y= 1/3 = const as a 
stiff-fluid solution. It is found that this solution is stable for early times (collapsing phase) but unstable for late times 
(expanding phase). This implies that there exists kink instability at the particle horizon scale for the flat Friedmann 
scalar-field universe and also for the flat Friedmann stiff-fluid universe. 

The Roberts solutions are exact and general self-similar solutions for the scalar-field system with k = 0, which will 
be described in Appendix According to the stability criterion obtained here, these solutions do not suffer from 
kink instability for early times while they do for late times. Since this solution has a null gradient of the scalar field 
at the similarity horizon, the stability criterion for self-similar stiff-fluid solutions is not applicable to this solution. 

We can find the static self-similar solution with a scalar field for Attk^ = 1, which is given by 

g = 2x, g = 4:X, ft. = ±(47r)"^/^ln|a;| -Fconst. (5.1) 
Also with a stiff fluid, we find the static self-similar solution as 

M=^ay^ S = al/\ V^laZ'^', 2/ = !, (5.2) 

where the regularity condition (|3.2U|) at the sonic point has been considered. Apparently, this solution seems unstable 
for early times but stable for late times. We can show, however, that this solution is always unstable. To demonstrate 
this explicitly, we change the coordinate system as 

w = -21n|M| +lni?, (5.3) 

for the scalar-field solution, and 

w = \n\t\, R^all^r, (5.4) 

for the stiff-fiuid solution. Then the line element for both cases is given in the static form as 

ds^ = -2R^du? + 2dR^ + R^{de + sin6|2rf0). (5.5) 

We should note that each of early-time and late-time solutions corresponds to the whole spacetime of the static solution 
through Eq. H5.3|l or (|5.4() . Inversely, when Eq. (|5.5() is given, it can be regarded as an early-time solution using the 
coordinate transformation Eq. H5.3() or (|5.4(l . In this solution, the similarity horizon condition g = 2x or sonic point 
condition = 1 is satisfied everywhere. It implies that we can insert a kink mode everywhere. Then the evolution 
equation (|2.3U|) or (|3.37|) for the kink mode perturbation holds along the radial null curve on which the discontinuity 
is inserted. Therefore, when we insert a kink mode perturbation, which vanishes inside and is discontinuous at some 
radius, it will blow up to infinity along the ingoing null curve {w + hiR = const) as time proceeds. In this sense, this 
solution is unstable. The situation is somewhat different for a perfect fluid with P = kp for < fe < 1, which will be 
described in Appendix lUl 

D. Horizons 

For self-similar stiff-fluid solutions, a sonic point may be a black hole event horizon, particle horizon or Cauchy 
horizon, simultaneously. Therefore, we can discuss the stability of these horizons. These horizons appear for late- 
time solutions. The black hole event horizon is characterized by an anti-transsonic point for late-time solutions. 
According to the result obtained here, ys > 1 must be satisfied on the event horizon. Then the black hole event 
horizon corresponds to a saddle and is stable against the kink mode. The particle horizon and Cauchy horizon are 
characterized by transsonic points for late-time solutions. On these horizons we can have < j/g < 1/2 or 1/2 < t/s < 1. 
These horizons correspond to nodes and are unstable if < j/s < 1/2 while they are stable against the kink mode if 
1/2 < Us < 1. The above also applies to similarity horizons in equivalent scalar- field solutions using correspondence 
relation H4.11|l . 
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E. Other applications 



Nonanalytic but regular self-similar stiff- fluid solutions are possible only for 1/2 < ?/s < 1. It is found that these 
solutions are unstable for early times but stable for late times. Anti-transsonic self-similar stiff-fluid solutions are 
possible only for 1 < y^. These solutions are unstable for early times but stable for late times. Again, the above also 
applies to equivalent scalar-field solutions using correspondence relation H4.11|l . 



VI. SUMMARY 



We have investigated the stability of self-similar solutions with a scalar field and those with a stiff fluid in general 
relativity. The kink instability, which we have considered here, was studied in a Newtonian gas system and in the 
Einstein-perfect-fluid system with the equation of state P — kp (paper I). Since only the fluid system was considered 
in these previous works, it was not clear whether the kink instability is unique to fluid dynamics or not. The present 
work, for the flrst time, has shown the existence of kink instability in self-similar solutions of the Einstein-Klein- 
Gordon system. The most intriguing feature is that the kink instability grows exponentially in terms of T = — ln|u| 
for the scalar-fleld system and also in terms of r = — In |t| for the stiff-fluid system, while it grows more rapidly and 
blows up to infinity at a finite moment before t = for the perfect-fluid system with P — kp [Q < k < 1). In other 
words, the kink instability results in the reduction of the rank of differentiability and the formation of a shock wave 
before the singularity formation for a perfect fluid with < fc < 1, while it does not for the scalar-fleld system and 
also for the stiff-fluid (fc = 1) system until the singularity forms. 

These systems have recently attracted attention in the context of critical phenomena in gravitational collapse. The 
present result shows that both the BCGN scalar-field solution and Evans-Coleman stiff-fiuid solution are unstable 
against the kink mode. This implies that the latter solution, which was identified with a critical solution in the 
stiff-fiuid collapse by Neilsen and Choptuik and Brady et al. cannot be a critical solution once we allow 
sufficiently small discontinuity in the density gradient field in the initial data sets; nor can the BCGN scalar-field 
solution be a critical solution once we allow sufficiently small discontinuity in the second-order derivative of the scalar 
field in the initial data sets. As another important application, we have shown the kink instability at the particle 
horizon of the fiat Friedmann universe with a scalar field and with a stiff fluid, while the flat Friedmann collapse 
solution is stable against this mode. 
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APPENDIX A: ROBERTS SOLUTION 

In the Bondi coordinates (|2.2|l . general solutions to Eqs. (|2.13|) - (|2.17|) with k = are given by |34 



(Al) 
(A2) 



2 



In ^. (A3) 



X 



where x is defined by Eq. (|2.8|l . and the solutions are parametrized by a and referred to as Roberts solutions. This 
solution has been discussed in the context of cosmic censorship [s^l and critical behaviour [3^ ^oj . 

For (7 = this solution reduces to the Minkowski spacetime and we restrict cr to cr ^ below. Then, these solutions 
admit no timelike Killing vector. From Eq. H4.2|l . the energy density of the equivalent stiff fluid is formally given by 



1 a2(2V^^T^-l) 
'-8^ x^ • 
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We can see that the solution has a spacetime singularity at a; = 0. If and only if < < 1/2 this solution has a 
region around the centre given by 



< a; < xi = - cr2^ (A5) 

where the gradient of the scalar field is spacelike. Ai x — xi the solution has a null gradient of the scalar field. If and 
only ifO<|(T|<l/2, the solution has a similarity horizon at 



(A6) 

Therefore, for < |<j| < 1/2, the hypersurface x = xi where the gradient of the scalar field is null coincides with 
the similarity horizon x — Xs- This implies that the stability analysis for self-similar stiff- fluid solutions developed in 
Sec, nil Cl cannot be used, although that for self-similar scalar-field solution developed in Sec. Ill Cl is still applicable. 
See for spherical and nonspherical perturbation with analyticity on the Roberts solution. 

APPENDIX B: HIGHER-ORDER DERIVATIVE TERM 

Let / be continuous on [0, oo) and continuously differentiable on (0,cxd) and have the limit value Ivcux^+oxf {x), 
where the prime denotes the differentiation with respect to the argument. When we define g by g{x) = xf{x), g is 
continuous on [0, oo) with g(0) = and continuously differentiable on (0, oo) and have the limit value lim^^+o g{h)/h = 
/(O). Then, from the mean value theorem, for any a > 0, there exists ^ (0 < ^ < a), such that 

(Bl) 

a 

In terms of /, we have 

f{a)=a'{0 + f{0- (B2) 

When we take the limit a +0, ^ also goes to zero. Since the left-hand side and the second term on the right-hand 
side approach the same limit value /(O), the first term on the right-hand side goes to zero in this limit. This implies 

lim xf'{x) = 0, (B3) 

because we have assumed the existence of the above limit. 

Applying the above proof, we can show that the right-hand side vanishes in Eq. (|2.29|) for the scalar-field system, 
where / corresponds to 6j', and in Eq. H3.35|l for the perfect- fiuid system, where / corresponds to Sr]'. 

APPENDIX C: STATIC SELF-SIMILAR PERFECT-FLUID SOLUTION 

We briefly review the static self-similar solution for a perfect fluid with the equation of state P — kp for < fc < 1. 
The case of fc = 1 is described in Sec. IV CI We use t and r as the coordinates given by Eq. ^i.'6\ . For < fc < 1, the 
line element in the static self-similar solution is written as 

ds^ = -alR^dw^ + PldR^ + R^{de^ + sin^ ed(j)), (Cl) 

after the following coordinate transformation 

{1 + fc . .iiiii 
-{-1)^+" +wq for t < 0, 
1 1 fc 1. (C2) 
-ti+fc +wq for t > 0, 
1 — fc 

R = ikv, (C3) 

where a/c, Pk and 7^ are some positive constants which can depend on Oq-, and fc, and wq is an arbitrary constant. 
In the above, w increases from —00 to +00 as t increases from —00 to +00. Therefore, the static self-similar solution 
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consists of both the early-time solution and the late-time solution. We can identify an arbitrary value of w with t = 
because of the arbitrary constant wq . Therefore any value of the time coordinate w can be regarded as in early time 
and also in late time. This is essential to the instability of the static self-similar solution as is discussed in paper I. 
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